Non-Markovian decay beyond the Fermi Golden Rule: Survival Collapse of the 
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The decay of a local spin excitation in an inhomogeneous spin chain is evaluated exactly: I) It 
starts quadratically up to a spreading time ts- II) It follows an exponential behavior governed by 
a self- consistent Fermi Colden Rule. Ill) At longer times, the exponential is overrun by an inverse 
power law describing return processes governed by quantum diffusion. At this last transition time tn 
a survival collapse becomes possible, bringing the polarization down by several orders of magnitude. 
We identify this strongly destructive interference as an antiresonance in the time domain. These 
' general phenomena are suitable for observation through an NMR experiment. 
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I. INTRODUCTION 



A typical quantum exponential decay involves a finite set of states in presence of an "environment" , i.e., weakly 
coupled to a set of states whose spectrum is dense. The decay of these states is usually described with the Fermi 
Golden Rule (FGR). However, this description contains approximations J3| that leave aside some intrinsically quantum 
' behaviors. Various works on models for nuclei, composite particles excited atoms in a free electromagnetic field 

J> I and in photonic lattices [1] , and models for decoherence , showed that the exponential decay has superimposed 
\^ . beats and does not hold for very short and very long times, compared with the lifetime of the system. The short time 
regime has received recent attention in connection to the Quantum Zeno Effect 0, [13, El, [13] and has been observed in 
,; trapped atoms [l^. In contrast, although different models predict some form of power law for long times [13. [isl [l6t . 
■ the cross over to this long time behavior has been neither experimentally observed nor physically interpreted. 

In this letter, we present a model describing the evolution of a local excitation in the otherwise homogeneous 
' polarization of a system of interacting spins. This situation has two desirable properties: 1) The full dynamics can 
be solved analytically and interpreted; 2) An actual Nuclear Magnetic Resonance (NMR) experiment can be tailored 
to observe this dynamics. Specifically, our model describes a linear chain of nuclear spins interacting under an XY 
(planar) interaction. In this situation, the evolution of a local spin excitation reduces to the dynamics of a localized 
density excitation in a system of non- interacting fermions [T^ . llSl . Il9l | . This excitation decays into a well resolved 
' wave packet propagating along the spin chain. Such decay could be observed with NMR because it also describes the 
^ ' dynamics of a multiple quantum coherence experiment [20| in a chain of spins with dipolar interactions in the solid 
O^' state ^3]. Furthermore, a full experimental dynamics of an effective XY Hamiltonian is achieved using NMR pulse 
sequences in liquid samples where the spin wave dynamics has been observed ,22,]. 

Based on our model, we are able to quantify the quantum nature of the deviations from the Fermi Golden Rule. We 
identify three well defined time regimes: 1) For short times the decay is quadratic {l—[tVo/h]^), as is expected when the 
coupling Vq of the local state with the continuum is treated perturbatively. This lasts for a time ts ~ hnNi (sr), where 
A^i (er) is the spectral density of the final states at the resonance energy e^; 2) An intermediate regime characterized 
by an exponential behavior, the self-consistent Fermi Golden Rule (SC-FGR) where the rate, the pre-exponential 
factor and the resonance energy are found self-consistently; 3) A long-time regime in which the exponential law is 
overrun by an inverse power law which is identified with the quantum diffusion in the chain. At this last cross-over, 
the oscillations could lead to a dip of several orders of magnitude in the local polarization. This survival collapse is 
identified with a destructive interference between the pure survival amplitude, i.e., the SC-FGR component, and the 
return amplitude, associated with higher orders in a perturbation theory. This striking quantum phenomenon can 
be seen as a dynamical version of the antiresonance that has been described for steady state observables [l^, [13, [1^ ■ 
Now, the destructive interference is also due to the splitting of the wave among two different families of pathways in 
space. However, the interference is now restricted to the narrow time window when the amplitudes are comparable 
and the phases are opposite. 
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II. DYNAMICS IN A NUCLEAR SPINS CHAIN 



We use the known mapping between spins and fermions [26[ together with a new formulation for spin dynamics 
based on the non-equihbrium Keldysh formahsm developed in Refs. [HI and [28| . 

The two spin correlation function in a system with M spins 1/2 evolving under a Hamiltonian H, 



Pf,^ (t) 





Sf it) Sf (to) 











(1) 



gives the amount of the z component of the local polarization on the site /th at time t, provided that the system 
was, at time to < t, in its equilibrium state with a spin 'up' added at the ith site. Here, (t) =e'^*5'|e~'^* is the 

spin operator in the Heisenberg representation and |^cq) — '^p^sn ^'oq' / is the many-body equilibrium mixed state 



oq 

constructed by adding states with different number N of spins up with the appropriate statistical weights and random 
phases. We will assume the high temperature limit, that leads to equal statistical weights |sAr|^ = ^{^)- 

We consider a linear chain of M spins in an external magnetic field. They interact with their nearest neighbors at 
distance a,, through an XY coupling: 



qy qy 



(2) 



where S'^ {u = x, y, z) represents the Cartesian spin operator. The first term of this Hamiltonian is the Zeeman energy, 
where ri„ is the chemical shift precession frequency; the second term, Hxy, contains the J„+i.„ coupling between 
sites n and n + 1. It gives the flip-flop interaction S^j^iS^ + S^j^^S^ in terms of the rising and lowering spin operator 

The Jordan- Wigner (J-W) transformation i2g| establishes the relation between spin and fermion operators at each 
site n. When H commutes with the number operator, the different subspaces N are decoupled. Further simplification 
is obtained for Hamiltonians which are quadratic in the fermionic operator, as the case of Hxy- Due to the short 
range interaction, after a J-W transformation, the only non-zero coupling terms are proportional to c^_|_]^c„ = 5'^^_]^5'~, 
where c^, c„ are the creation and destruction operators for fermions. The Hamiltonian become 
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where e„ = hQ.n are the site energies and Vn 



+ l,n 



\Jn+i.n are the hoppings. Each subspace has N non-interacting 



fermions. The eigenfunctions ^'^^^^ are expressed as a single Slater determinant built-up upon the single particle 

wave functions "0^ describing a particle of energy in a chain. Under this condition and defining \i) = |0) , with 
|0) the fermion vacuum, Eq.([T]) reduces to 

Pf,r it) 
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where (t) is the retarded Green's function for a single fermion. 

Therefore, for systems represented by a 1-d chain of spins with nearest neighbors XY interaction, at high tem- 
perature, the dynamics of a local polarization amplitude corresponds exactly to the wave function of single particle 
evolving according to a tight-binding Hamiltonian. 



III. LOCAL EXCITATION: THE EXPONENTIAL DECAY AND BEYOND 

Let us describe the evolution of a local excitation |0) = Cfj" |0) in a Hamiltonian whose spectrum has a finite support. 
This is the case of most excitations in a lattice. The autocorrelation function is Eq.(I3]) with \i) = |/) ~ |0) . Expanding 
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the initial condition in the eigenstates \ipk} one obtains [1,0] 
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The term in brackets is the Local Density of States (LDoS) A^'o (e) at site 0th. It can be evaluated using the retarded 
Green's function, 
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iVo(£) =--Im / dt G«(i)e"* 
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Then, we can express the autocorrelation function as the Fourier transform of the LDoS, 



Poo{t) - 
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This expression has numerical and analytical advantages because the Green's function can be accurately calculated 
in the energy representation, and the integral is limited to the spectral support. Besides, a clear identification of 
quantum interferences will be obtained by analyzing the argument under the modulus operator. 
Alternatively, the autocorrelation function can be written as 



Poo{t) = 0{t) / duj Muj) exphiwt], 

^ — OO 

where the spectral density of the particle excitations at site 0th: 

poo 

Jo{uj) = h de NQ{e)No{e + hoj), 
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has a direct physical interpretation and can be easily computed [29|. 

All the previous equations remain valid when the size of the system, and hence the dimension of the Hilbert space, 
becomes unbounded {M — > oo). In this case, either part or the whole of the discrete (pure point) spectrum, may 
become a continuous energy band of delocalized (extended) states in the finite range [eL,£c/]- If the system does not 
present localized eigenstates ^3], ^ois) vanishes outside the band (Fig. [T|). On the other hand, if the initial state |0) 
has a finite weight over one or more localized states its evolution can not fully decay. Here, we consider cases that 
exclude such situation. Hence, if |0) requires an expansion in an infinite number of eigenstates, its evolution becomes 

an irreversible decay. In particular, the unperturbed state of energy eq — (o\H\0^ becomes a well defined resonance 

if |0) can be expanded in terms of the eigenstates within a small breath Fq around an energy = + ^Oi where Aq 
is a small shift due to the interaction. Furthermore, the validity of the Fermi Golden Rule for Poo(0 requires fsTl that 
the state |0) is similarly coupled to each of the unperturbed states |^^) with energies in a continuum spectrum. 

In order to evaluate the local dynamics, we perform the integral in Eq.([7]) using the residue theorem and following 
the path shown in the Fig. [TJ In the analytical continuation No{z) = No{e + ie'), resonances appear like poles in 
the complex plane. We will consider Hamiltonians where an initially localized state with energy eg interacting with 
a continuum gives rise to a single resonance, i.e., the LDoS presents poles at Sr ± iFo. The van Hove singularities on 
the contour are excluded with circle arcs with radii R. Their contribution to the integral vanish when i? — > 0, because 
the band edges are of the form (e — e^Y with v > —1. Also, the integral over the contour z — e — iL; e e [eL,Su], 
vanish when L oo. Then, we obtain 



Poo(t) = [a e-(^^+'"'-)*/^, + / de'e-^'*/'^[e-'^^*/'"Wo(£L - ie') - e-'^^*/'*iVo(£c/ - ie')]!', (10) 

SC-FGR 



return correction from quantum diffusion 
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FIG. 1: Local spectrum (LDoS) in the complex plane z — e+i e. el and eu are the lower and upper band-edges, respectively. 
£i and £2 are localized states. The resonance energy is Sr ~ eo — Aq, and the pole appears in Sr — iFo. The integration path is 
shown with dotted lines; consist of four straight lines and two arcs, that avoid the band-edges singularities. 



where a — limz^e^,_iro [27ri (z — + iFo) Nq{z)] and t > 0. If we approximate the LDoS by a Lorentzian function 
that jumps to zero outside the band, we can see that 



A=\af~l + S, (11) 



where 



0<(5=-- < 1. (12) 

TT (e,. - El) [eu - £r) 

The first term of Eq. (jlOp already supersedes the usual Fermi Golden Rule approximation since it has a pre-exponential 
factor [A > 1) and the exact rate of decay Fq. This result is the self- consistent Fermi Golden Rule (SC-FGR). By 
analogy with a classical Markov chain, this exponential term is identified with a "pure survival" amplitude. Within 
the same analogy, the second term will be called "return" amplitude, as it is fed upon the initial decay. The first 
is the dominant one for a wide range of times, while the diffusive decay of the second, dominates for long times 
and brings out the details of the spectral structure of the system. In the quantum case, the second term is also 
fundamental for the normalization at very short times where the most excited energy states of the whole system can 
be virtually explored. Both terms combine to provide the initial quadratic decay (Quantum Zeno regime) required 
by the perturbation theory: 

^oo(0-l-^((£-£.)')^„ + --- , (13) 
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-l4^-')^o + ---- (14) 

Here ((e — e^)^)^ and {t^'^) j are the energy and frequency second moments of the densities Nq (e) and Jo{u!), 
respectively. This expansion holds for a time shorter than the spreading time ts of the wave packet formed by the 
decay. In other systems, the divergence of the second moment leads to different short time decays ^32]. 

For long times, the behavior of Pqo (t) is governed by the slowly decaying second term in Eq. fTU)) . Only small values 
of e' contribute to the integral. In turn this restricts the integration of the LDoS to a range near the band-edges. 
Then, one can go back to Eq.Q and retain only the van Hove singularities of the Local Density of States at these 
edges (e.g. No{e) oc 9{e — el) (e — £lY which implies [1^ j7o(w) oc 9{uj)uj'^'^^^). Each singularity would contribute to 
the slow decay at long times (Poo(0 °^ jtP^'''^^^). The relative participation of the energy states at each edge of the 
LDoS is given by the relative weight of the Lorentzian tails at these edges: 

3- ier-e,f+Tl 
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Then, the polarization for long times is 



Poo{t) - [l + - 2l3cos{Bt/h)] 



de'e-^'*/'"Wo(eL -ie') 
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(16) 



where B = ejj — e^. This means that the long time behavior is just the power law decay of the integral multiplied by 
a factor having a modulation with frequency B /h. 

IV. SURVIVAL COLLAPSE 

In steady state transport [1^ as well as in dynamical electron transfer [2^ there are situations in which a particle 
can reach the final state following two alternative pathways. Since each of them collects a different phase, this allows a 
destructive interference blocking the final state. This phenomenon has been dubbed antiresonance [IS, [la]- It extends 
the Fano resonances, which describes the anomalous ionization cross-section [23|. In the present case, the survival of 
the local excitation also recognizes two alternative pathways: the pure survival amplitude which is typically described 
by the Fermi Golden Rule, and the paths where the excitation has decayed, explored the environment, and then 
returns. These two alternatives can interfere. We rewrite Eq.(fTO|) to emphasize that the local polarization Pqo [t) is 
the result of two different contributions: 

Poo {t) = l^-s + = |*s|' + + 2Re[^'J*fl], (17) 
where the phase in n arise from the exponentials with el and eu (the LDoS is real for any argument). Hence, 

(t) = \a\ e-i'^'-e-^"*/ V^^--^^)*/'', (18) 

*fl,(i) = l*«W|e''^^*^; (19) 

f P sin (Bt/h) \ 

(h (t) = arctan — \' ' . 20) 

\1- 13 cos (Bt/n) J ^ ' 

where Eg. ipO)) results using the long time limit of Eq. p^ . 

While the interference term in Pqq (t) is present along the whole exponential regime, it becomes important when 
both, the pure survival amplitude and the return contribution, are of the same order. This occurs at the cross-over 
time tn between the exponential regime and the power law. The interference term can produce a survival collapse, 
i.e., a pronounced dip that takes Poo (t) close to zero (see Fig. [2]). In order to obtain a full collapse, two simultaneous 
conditions are needed; 

|*5(ifl)| = |*fl(ii?)| and (21) 
{£r - eL)tR/h~ (j> {tn) = (vr - (/)a) + 27rn, n integer, (22) 

which are satisfied with a fair precision because 

|(£^ - el) /h\ > To/h > 2Tr /tR > \(t) {tR)\ /tR, (23) 

i.e., while the return amplitude has a phase with a slow variation, the pure survival term oscillates rapidly. When 
both amplitudes are of the same order, the destructive interference will be noticeable. 



V. DECAY IN A SEMI-INFINITE CHAIN 



Now we focus on a specific case of Eq.© that can be achieved experimentally and has simple analytical properties. 
We consider the Hamiltonian of Eq. ([3]) with the 0th site (spin) in the chain different from the others sites in both site 
energy (chemical shift) and hopping (J-couphng), i.e., eq ^ En = 2V and Vo,i = Vq < Vn,n+i = V for n > 0. This 
defines a continuous spectrum in the range [0, B = AV] which, in the lower edge, describes a particle of mass m in 
the continuum with V = /{2ma'^). Our model presents a resonance provided that the site ener gy i s not to close to 
the band edge, i.e., |eo — 2V\ < 2V — Vq /V. Otherwise, |0) would give rise to a locahzed state jSj]. The LDoS for 
this problem is evaluated using the Dyson equation 

\vo.iGf[,iE)Vi^o, (24) 
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following the general continued fraction procedure described in Ref. [33l |: 

1 « 1^1 H (41/ — ( -^'^ 

No{e 
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(25) 



Note that, because of surface effects in the seminfinite d-dimensional space, the LDoS has van Hove singularities of 

the form N['^\e) cx which differ from those in the bulk N'^'^^e) cx £(''-2)/2. pig. [T] shows 7Vo(e) for Vo/V = 0.4 
and Eq/V = 1. The resonant state (the poles of the LDoS) appears in Er ± iro, where 



Sr = So + Ao; Ao = 



1^2 
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(26) 
(27) 



Identifying the local density of states at the first site in absence of interactions with the 0-site as 



7Vi(e) = -ilmGi,i(e) 



16 1 , , 



(28) 



Note that OoF (e) /h is the group velocity of a wave packet with energy e and F (eq) ~ Fc. One realizes that the 
expression of Eq. ([25]) factorizes as a pure Lorentzian and -/Vi(e) : 



No(s) 



Fc (e 



F2 
^ 



(29) 



Then, applying the convolution theorem to Eq. ([7|) we get a convolution integral of two functions in the time domain 
with well characterized time dependence. 



(30) 



The first factor inside the integral is the renormalized survival amplitude as described by the SC-FGR. The second 
factor is the return amplitude to site 1 in a semi-infinite chain where site is missing. It is expressed in term of 
the Bessel function of the first kind as g (t) = 2e~^^^*/'^Ji{2Vt/h)/{2Vt/h), which shows fast oscillations and decays 
with the power law t~^^'^. This describes the quantum diffusion in the chain (28l. [33|. It appears convoluted with an 
exponential kernel whose oscillation and decay have a longer time scale. For positive times g (t) coincides with the 
response function. This knowledge allows us to solve the integral in the different time regimes (short, exponential and 
long time). After some algebra we get 



1 - (Vot/hY 



t < ts 



Pooit) 



C 



^exp(-2Foi/ft), ts<t<tR 



2/3 



sin (Bt/h) 



r(e.)t 



(31) 



,tR<t 



where ts is the cross-over time from the short time regime to the exponential SC-FGR, and time tj^ separates the 
SC-FGR and the power law regime. Also, 



c 
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F2 > 1 
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(32) 
(33) 
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0{V^/V) in Eq.dlll) to obtain Tq 



and (3 was defined in Eq. ((T5|) . Here, we used Vb <C F and eq ~ e 
which coincides with the SC-FGR. Near the band edge Fq ~ 9>Vq / B^Jsr/ B, and T{er) 
and averaging in a period, one gets 



Pm it) 



An Er 



Tier)t 



For long times. 



(34) 



At long times, the probability of finding the particle at site is proportional to the probability of finding it at site 1, 
i.e., Pooit) ~ (Vo/so)^ PQi{t). Hence, the factor gives the probability of tunnehng back to |0) . It meassures how the 
component of the band edge (that determine the long time behavior) over the surface state |1) mixes with state |0). 
The assignation of a time scale to the return probability in the last term is arbitrary. We choose F (e^) , the dominant 
group velocity of propagating wave packet of energy e^- Hence, the second factor becomes the inverse of the number 
of cycles within the main decay. 

It is important to note that the cubic power law decay obtained for long times is a consequence of the ^/e dependence 
of the band-edges of the LDoS, i.e., this power law is consistent with Eg. p^ taken together with Eq. pSj) . Notice 
also that the short time scale, H/Vq, can also be obtained from the local second moment of the Hamiltonian. 

From the analytical result given in Eg. lpTj) . we get the characteristic times ts and t^. A good estimate of ts is 
obtained from the minimal distance between the short time decay and the exponential: 



d_ 



- [Vot/hy - ^exp(-2Foi/ft) 
Expanding the exponential in its Taylor series, we get 



ts 



HTqA 



0. 



for Er < 2V. 



(35) 



(36) 



Here, we see that in this parametric regime the short time cross-over is only determined by Ni (e^) , the local density 
of states at the first site of the chain. We may invoke the optical theorem [s^, to interpret A^i (e^) as the time scale 
at which a wave packet with energy Sr esca pes from the 1st site region, i.e., the excitation build from the decay, 
preventing the return to the original 0th site [33l [35| . 

The time tji is obtained from the cross-over between the exponential regime and the power law decay 



^exp(-2FoW;i) = C 



T{er)tR 



(37) 



We can use yA/C ~ 2v^£o/Vb-\/e77rO' and solve iteratively the transcendental equation, i.e., 

,(0) 



t 



(n+l) 



Fo V ^ VnVTj 2Fo 



r(er) (n) 



(38) 



Already in the third order we get a very good agreement with the cross-over observed in the exact dynamics. 



VI. NUMERICAL VERIFICATION 



We verify the above results following two independent procedures. Since one has a closed analytical expression for 
iVo(£), the numerical Fourier transform is straightforward. Alternatively, we find the dynamics from the numerical 
eigenvalues and eigenvectors of the finite system with M sites. Both of them coincide as long as M is big enough so that 
the finite system effects become negligible. This requires that the mesoscopic echo [l^, arising at a time tME ~ hM/B, 
appears well beyond the cross-over time tj^. Both procedures provide perfect agreement with the analytical results. 
In Fig. [5] we show Pooit) in a semilogarithmic scale. The exact decay confirms the time dependences exhibited by 
the analytical approximation of Eq. (j3ip . The initial quadratic decay is amplified in the upper inset. Then, the curve 
is followed by the exponential SC-FGR. Finally, it presents a cross-over at to the asymptotic power law decay. 
This time-scale is easily identified through the survival collapse shown as a dip in the survival probability. There, 
the polarization suddenly decreases from its average by almost three orders of magnitude. The inset on the bottom 
shows the small oscillation that modulates the power law. 
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FIG. 2: Local polarization, in a semilogarithmic scale, as a function of time. We consider an unperturbed energy of EqIV — 1 
and interaction strength Vo/V — 0.4, that leads to a resonance energy Sr/V = 0.9 and an exponential rate Fo/V — 0.14. This 
is the case that we consider in Fig[T] The decay exhibits: 1) The quadratic perturbative regime, which is shown amplified in 
the upper inset. 2) The exponential behavior as described by the self-consistent Fermi Golden Rule. 3) An asymptotic cubic 
power law decay, where 6 = C/F (er)^ (Ea l33|l . The lower inset shows the oscillation that modulates this decay. The cross-over 
time t[i when the survival collapse takes place is indicated. 




FIG. 3: Local polarization, in a semilogarithmic scale, as a function of time, in units of fi/V, for eo/V = 1.3, Vo/V — 0.75, that 
leads to a resonance energy Er/V = 0.85 and a SC-FGR exponential, shown with a dotted line, with rate To/V = 0.72 and 
pre-exponential factor A = 2.86. The "environment" has M — 1 — 19 spins. The decay exhibits a noticeable survival collapse 
followed by a cubic power law modulated with a well defined frecuency. At later times, the mesoscopic echo shows up. 



Since the model solved in the previous section could be applied to spins in a molecule or excitations in a designed 
nanostructurc, both of which have finite size, it is interesting to verify that the main features discussed also could be 
observed in such situations. In Fig. Owe show the dynamics of one spin in presence of an "environment" consisting 
on a chain of 19 identical 1/2 spins. The three regimes just discussed are clearly manifested. Later on, it appears 
a mesoscopic echo at tV/h > 20. Note that already at tj^ « 6h/V the magnetization decreases in seven orders 
of magnitude. For a brief period around tn coherent interference ensure an almost complete depolarization of the 
surface site that could not be achieved through decoherent decay. 
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VII. CONCLUSIONS 



In the present work we have discussed the exact dynamics of a local excitation that decays through the interaction 
with a continuum spectrum with finite support that acts as an "environment" . Our approach goes beyond the usual 
Markovian approximation that uses the Fermi Golden Rule to describe these environmental interactions. Within a 
simple, yet realistic model of a linear chain of nuclear spins with XY interaction, we found the exact behavior of the 
autocorrelation function for all times. The evolution starts with the expected quadratic decay. Then, it follows the 
usual exponential FGR regime, but with a corrected rate and a pre-exponential factor, i.e., the SC-FGR. Finally, 
we get the long time regime, that consists of a cubic power law decay modulated by oscillations whose frequency is 
determined by the bandwidth. This power law decay is a consequence of the iVe behavior of the LDoS in the band-edge 
fEg. p^ ). A similar result is obtained in models for unstable nuclei [1, 0, @, an atomic excitation in the free space 
and interacting with a photonic band Q. In those cases, the decay law is the regular Van Hove singularity of 
the free space. Here, the surface modifies the singularity and hence, the time decay. Also, we found the analytical 
expressions for the cross-over times ts and Ir of Eq. (|3ip enabling us to assert the range of validity of each regime. 

Finally, we find and quantify the survival collapse. This effect, hinted but not explained in previous works, is 
visualized as the destructive interference between the pure survival amplitude and the return amplitude that arises 
from pathways that have already explored the rest of the system. This non-Markovian result fully considers the 
memory effects to infinite order. 

In summary, through the exhaustive solution of a particular model, we made a conceptual analysis of a general 
quantum decay process applicable to the great variety of systems where a quantum exponential decay is observed. 
Besides this generality, what gives a particular interest to our model is its suitabilit y fo r an experimental test. This 
would imply the same procedure devised [22] to test the mesoscopic echoes [13, llSj j . In order to tailor an XY 
Hamiltonian in an NMR experiment, it uses a radio frequency pulse sequence that produces the truncation of the 
natural Heisenberg (J-coupling) Hamiltonian. Alternatively, using the relationship between the dynamics described 



by an XY Hamiltonian and multiple quantum dynamics 21|, the dynamics of our model could be observed with 



multiple-quantum experiments in solid state NMR [20| . The application of one of the above procedures to relatively 
small linear molecules would enable the observation of the survival collapse. One could freeze the dynamics at this 
time obtaining an almost null survival of the local excitation. Since the survival collapse depends critically on the 
cooperative coherence of the whole system, it would be quite sensitive to decoherent processes and hence it could be 
applied to evaluate them. 
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